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ABSTRACT 

The CP violation effects in long baseline neutrino oscillations are studied in the framework 
of a four-neutrino model (three active neutrinos and one sterile neutrino). It is assumed that 
neutrino masses are devided into two nearly degenerate pairs, as indicated by the oscillation 
data. Approximative analytic expressions are derived for the probability differences AP 0( g = 
P{y a v p) — P(V a — > Vp) taking into account the CP violation effect of the Earth's crust. 
The matter effect is found to be small compared with the genuine CP violation term, in 
contrast with the three-neutrino model. 
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A new stage of the neutrino oscillation studies is represented by the long baseline (LBL) 
oscillation experiments. The LBL accelerator experiment K2K fl| begins to take data in this 
year (1999), whereas the MINOS || and a CERN-to-Gran Sasso project || will start in the 
first year of the next millenium. Some authors @, [|, |6|, [7], |8| have discussed the possibility 
of observing the leptonic CP violation by measuring in LBL experiments the difference of 
transition probabilities between CP-conjugate neutrino oscillation channels || pL0f| 3 such as 
Vn —> v e and V ^ — > V e . The direct measurement is, however, very difficult in the planned LBL 
experiments since the magnitude of the probability difference AP^ = P[y a —> Vp)—P{p a — > 
Vp) is usually expected to be below 0.01 and the difference of energy distributions of neutrino 
beams and disturbs this measurement at the level of (9(0.01). Moreover, the matter 
effects due to the Earth's crust, while not that important for the oscillation probabilities 
themselves, can be sizeable for the probability differencies and make it difficult to extract 
from the data the genuine CP violating effect caused by complex elements of the neutrino 
mixing matrix. 

On the other hand, a neutrino factory, a high intensity and high quality neutrino beam 
from a very intense muon source, has recently been under discussion in connection with the 
muon collider studies [O], 0. It would provide an excellent possibility to search for CP 
violation in LBL neutrino oscillations because the flavour content of the beam would be 
known (50% and 50% V e in the case of negative muons), as would be also the energy 



distributions of the neutrinos as soon as the muon energy and polarization are given [|12| . 

The magnitude of the CP violation in the neutrino oscillation depends on the number of 
neutrinos involved. One of the present authors (M.T.) has studied prospects for searching 
the CP violation at a neutrino factory in the scheme of three ordinary neutrinos z/ e , 
and u T , and he proposed a method for extracting the genuine CP violation effect from 
the observed signal [pT3| . It is not guaranteed, however, that the three left-handed active 
neutrinos are the only neutrinos that exist in nature. Sterile neutrinos, neutral leptons that 



2 



lack the Standard Model interactions, may also exist and they may give rise to interesting 
phenomena [|TJ], [TJ|, [16], |17], |18| . Actually, in order to explain the solar [|19j , atmospheric 



and LSND |2T| observations in terms of neutrino oscillations one necessarily needs at least 
one sterile neutrino in addition to the three ordinary neutrino flavours. 

In this paper we shall study the CP violation effects in the oscillations of four neutrinos 
and their signatures in LBL oscillation experiments at a neutrino factory. We assume that, 
besides the three ordinary neutrinos z/ e , and v T , there exists one sterile neutrino v s . We 
shall make a comparison of the predictions for the CP violation in the three-neutrino and 
four-neutrino cases and discuss how these two cases could be distinguished experimentally. 



In order to be compatible with the solar [|T^], atmospheric P0| , LSND data |j21| , as well 
as with the other accelerator and reactor limits, the neutrino masses in the four-neutrino 
model should be divided into two pairs of nearly degenerate masses [17 , |18|| , e.g., Am 2 & = 
Am^ <C Am^ tm = Am^ <C Am 2 LSND = Am^, where Am^ = m 2 — m 2 and m ; is the mass 
of the massive neutrino state [i = 0, 1,2,3). One can assume without loss of generality 
that < mo, mi < rri2 < m%. If the solar neutrino oscillations are driven by the MSW effect 



22|, one must require itlq > mi. In the following we will assume the four- neutrino spectra 



to be like this. 

Neutrino oscillation probabilities P{v a — > are determined by mass-squared differences 
Amjj and the elements of the neutrino mixing matrix U that connects the massive neutrino 
states z/j and the flavour neutrino states v a through the relation v a = Y^iVai^i p3| - The 
most general mixing matrix U for four Majorana neutrinos can be parametrized in terms of 
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6 rotation angles and 6 phases as follows: 
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where c y - = cos^ y - and s y - = smO^e 1 ^ |24" |. If the mass eigenstate neutrinos are Dirac 
particles, only three phases are physically meaningful. 

The probability of the neutrino flavor oscillation v a — > (a, /3 = s, e, fi, r) in vacuum is 
given by 

P(v a -> u 
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Am| = m 2 k — m 2 



L is the oscillation distance, and E is the neutrino energy. The probability P(v c 
the CP conjugated channel is obtained by replacing U with U*. 

The difference of the transition probabilities between CP-conjugate channels, 

AP af3 = P{v a -> Z/g) - P(17 Q -> 77^) 

= 4(J°J sin 2A 10 + sin 2A 20 + sin 2A 30 
+ sin 2A 2 i + sin 2A 3i + J 3 sin 2A 32) > 



(2) 

(3) 
(4) 
Vp) of 



(5) 
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directly measures the CP violation that originates from the complex phases of the mixing 
matrix U || [Hj. Basically due to unitarity of the mixing matrix U, only three of the dif- 
ferences AP a p are independent in the four-neutrino case, and, correspondingly, only three 
out of the six phases of the mixing matrix U of Majorana neutrinos can be determined by 
neutrino oscillation measurements. Therefore, as far as neutrino oscillations are concerned, 
our considerations apply to both Dirac and Majorana neutrinos. The three remaining in- 
dependent phases of U enter into the mass matrix elements and appear in lepton number 
violating processes involving Majorana neutrinos, such as neutrinoless double beta decay. 
By using the relations 

t12 _ t23 _ t02 t13 _ _ t23 _ t03 / r \ 

J a(3 — J a/3 J a/3 > J a/3 ~ J a/3 J a/3 > \ u ) 

one can rewrite (|5]) in the form 

AP a/3 = -4 Jg(sin 2Ai 2 + sin 2A 23 + sin 2A 3i ) + 4 sin 2 Ai 

+ 4J°J (sin 2A 20 - sin 2A 21 ) + 4J$(sin 2A 30 - sin 2A 3 i). (7) 

Since 

sin 2A 2 o - sin 2A 2i = 2 cos(A 20 + A 2i ) sin A 10 , 

sin 2A 30 - sin 2A 31 = 2 cos(A 30 + A 31 ) sin A 10 , (8) 

one can see that the last two terms in eq. (0) average out to zero in the LBL experiments 
due to the rapidly oscillating cosines. 

In the neutrino factory one has a + V e (z7 M + i/ e ) beam originating in the decay of 
fi~ As mentioned above, there are only three independent probability differences 

AP a p from which information about the CP violation can be inferred. Let us consider the 
probability differences associated with the channels v e — > v^, — > i/ e , and u e — > u T , which 
are observable in practice. They are given as 

AP efl = — 4J^(sin 2A 12 + sin 2A 23 + sin 2A 31 ) + 4 sin2Ai , 

5 



AP Me = -4jj 3 (sin 2A 12 + sin 2A 23 + sin 2A 3 i) + 4 J^ e sin 2A 10 , 

AP eT = -4 J 23 (sin 2A 12 + sin 2A 23 + sin 2A 3 i) + 4j£ sin 2A 10 . (9) 

In order to be able to estimate J%, we need information of mixings between the ster- 
ile neutrino u s and left-handed active neutrinos u e , z/ M and v T . These mixings are severely 
constrained by the standard Big Bang cosmology. If the mixing is too large, neutrino oscilla- 
tions, acting as an effective interaction, would bring the sterile neutrino in equilibrium before 
neutrino decoupling, and the resulting excess in energy density would endanger the standard 
scheme for the nucleosynthesis of light elements (BBN) p5| . Moreover, a new mechanism of 
resonant sterile neutrino conversion in the boundaries of spatial opposite-sign lepton number 
domains was recently proposed |26| and was found to lead to even more stringent constraints 
on the mixings, in particular for large mass-squared differencies. The ensuing constraints 
are summarized as follows P5| : for the v^ T — v a mixings, 

|Am 2 |sin 2 2# < 7 x 10~ 5 eV 2 for |Am 2 | < 2.5 x 10 3 eV 2 , 

sin 2 29 < 3 x KT 8 for |Am 2 | > 2.5 x 10 3 eV 2 , (10) 

and for the v e — v s mixing, 

|Am 2 |sin 2 2fl < 5 x l(T 8 eV 2 for |Am 2 | < 4eV 2 , 

sin 2 29 < 1(T 8 for | Am 2 | > 4 eV 2 . (11) 

From this we conclude that the z/^ jT — z/ s mixings can be neglected in the following study as the 
corresponding mass-squared differences are supposed to be of the order of Am1 SND ~ 1 eV 2 . 
On the other hand, for the v e — v s mixing, which is assumed to be responsible for the solar 
neutrino deficit, the mass-squared difference is very small, so that ([TT| ) does not constrain 
that mixing severely. 
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Neglecting the v^ T — v s mixings, the matrix U is expressed as follows: 

/ Q)l Sq^ S n 9 Sna \ 
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By using this approximation, we have 



J 2 1 = -Jll = - J 23 . = C 23 Sl2S23Sl3Sin0 = J CP , 



(13) 



where (ft = 0i 3 — <pi2 — 023- Jcp is the four-neutrino counterpart of the famous Jarlskog 
invariant E7[ defined for the three flavour mixing. Then, eqs. (0) become 



efM = -AJ CP (sm 2A 12 + sin 2A 23 + sin 2A 31 ) + 4j£ sin 2A 10 , 
\ e = 4 J CP (sin 2A 12 + sin 2A 23 + sin 2 A 3 i) + Uf e sin 2A 10 , 
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4 J CP (sin 2A 12 + sin 2A 23 + sin 2A 31 ) + 4 j£ sin 2A 10 . 



(14) 



Taking account of = — J®1, we get the CPT invariant relation AP eAt = — AP^ e , which is 
the same relation as in the three family model. 

How large is J^? In the approximate mixing matrix U, that is when sq2 = and S03 = 0, 
we have = = 0. Therefore, the CP violating effect follows from the active neutrino 
mixings and phases. The mixing between the sterile and the active neutrinos contributes to 
AP a p only indirectly through the unitarity of the 4x4 mixing matrix. 

The above-derived CP violation measures AP a/ 3 as such should not be compared with 
observations but one should take into account also matter effects in the earth. The earth is 
CP-odd in the sense that it acts differently in the propagation of neutrinos and antineutrinos. 
It will therefore affect AP a p, even if the distance travelled by neutrinos is less than 1000 Km, 



like it is in LBL experiments. The matter effect should be carefully analyzed since it depends 
strongly on the mass hierarchy, mixings and the incident energy of the neutrino, as was shown 
in the previous works |28|, [29], |3(J. The effective mass squared in the matter, M^, for the 
neutrino energy E is in weak basis given by 
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V / 2Gi?n e /2. For antineutrinos, the effective mass squared is 
a' and U -> U*. 



—a, a - 

Since Am\ <C a, a' in LBL experiments, the sterile neutrino decouples from active ones 
as far as S02 — S03 — 0. Therefore it is enough to discuss the matter effect on the active 
neutrinos only in LBL experiments Q 

Although we can calculate AP a p in the presence of matter numerically, an approximate 
analytic formula is useful for investigating the qualitative structure of the matter effect. For 
the case that the highest squared-mass difference is 0(leV 2 ) the appropriate formulae have 
been given in the lowest order approximation by Minakata and Nunokawa M: 
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4 In short baseline experiments and the solar neutrino experiment, the sterile neutrino couples with active 
ones as far as sqi ^ 0. 
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and SV denotes the first order correction for the mixing matrix in vacuum as 
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The last term in ([If]) is at least second order in a. The quantity Iij(a) appearing in fllq) is 
given by 
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Eq. (|16D relies only on the hierarchy a <C \Amfj\/E but is not affected by any possible 
hierarchy among Am 2 j. 

We now approximate the density of the Earth by a constant (p = 2.8 g/cm 3 ), i.e. we take 
a as a constant fj By using eqs. fllTf ) and (0), we get | 
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(20) 



where the last term of eq. ([lq ) proportional to lm(UUU5V) has been ignored because it is 



of order a 2 or higher. The AM 2 and Am 2 denote Am 2 x ~ Am 2 x ~ (9(eV 2 ) and Am 



32 



10 3 eV 2 , respectively. Since AM 2 ^> Am 2 , the terms proportional to a/ AM 2 can be safely 
neglected in (^0]) . The first term in (|20| ) is the genuine CP violation and the second and third 



terms are leading matter corrections. We can write approximate formulae for the u e 



v. 



5 The constant matter density is not always a particularly good approximation quantitatively, so that 
for an accurate estimation of the matter effect one should apply a real earth model with an appropriately 
varying density. Actually, Koike and Sato have discussed the matter effect in the K2K experiment by using 
a real earth model |H| . 

6 The case of Am^ ~ Am| 2 ^> Amf 2 was discussed in rcf.||. 
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v e — > v T and — > z/ e channels in terms of the mixing parameters as follows: 
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Taking account of s 12 , s± 3 <^ 1, these formulae are expressed as 
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We have checked that the approximate results in eq. (24) are in agreement with numer- 
ically computed results within 1% for the parameter values L = 732 km and E = 7GeV. 
According to eq. fl25|) the matter effect is of the order of s 2 2 s 2 3. The mixing parameters S12 
and S13 are constrained by the accelerator and reactor data. An upper limit s 2 2 .s 2 3 ^ 10~ 4 
can be inferred from the results of the Bugey disappearance experiment [0 . By inserting 
s 23 ~ I/a/2, as indicated by the atmospheric neutrino data, one obtains \Jcp\ ^ 0.005 sin (f> 
and thus, by using the relevant parameter values L = 732 km, E = 7 GeV and p = 2.8 g/ cm 3 , 
one has \K matt \ < 4x 1CT 5 and \ J ma tt\ ^ 10~ 5 . On the other hand, the genuine CP violation 
4J C p sin(Am 2 L/2P) could be as large as 0(iO -2 ). Thus the matter effect is quite tiny in the 
four-neutrino model, in contrast with the three-neutrino model where it may with certain 



parameter values be as important as the genuine CP violation term (13 
From (El) one obtains the relations 
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(26) 



The first relation shows that in the four-neutrino model the vacuum relation AP efl = —AP eT 
is not spoiled by matter effect, in contrast with the three- neutrino case where in general 
|AP er | > |AP eM | and for some part of the allowed parameter space even |AP er | ^> |AP e/i | 
due to the matter effects (if we assume the sine of the phase angle to be positive; in the 
opposite case the relations turn around Jl3|). The matter independence of the second relation 



is guaranteed by CPT invariance. In the three-neutrino model [13] one has 



AP eM ~ W CP fcP + P m s 
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where J' CP is the Jarlskog parameter in the three-neutrino case, 

o / 4aP 
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Am?,L Amli AmLL , , 

/cp = sm _^ + sill _^ + sm _|_, (28) 



and two independent mass differences are —Am\ 2 = Arrig ~ 10 5 ~ 10 4 eV 2 and Am| 
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Am^ m ~ 10~ 3 eV . In formulas fl27"D the value of J' CP , unlike Jcp in the four-neutrino 
case, depends on the solution of the solar neutrino problem. Because of the differences in 
equations (|24]) and (^7|), it may be possible to distinguish the three- neutrino model from the 
four-neutrino model by measuring the quantity AP a p in all three channels. 

In summary, we have studied the CP violation in four-neutrino oscillations in a long 
baseline at a neutrino factory. We have evaluated the probability differences AP a p = 
P{y a — vp) — P{v a — Vp) between CP-conjugate channels and compared them with the corre- 
sponding results in the three-neutrino case. For the analytical expressions of these quantities 
we have used an approximation in which the mixing of the sterile neutrino with the muon 
and tau neutrinos are neglected. This approximation is justified by cosmological constraints 
on the active-sterile mixings. The matter effect on the CP violation has been derived fol- 
lowing the approximative approach of Minakata and Nunokawa ||. We conclude that in 
the four-neutrino model the matter effect is generally small compared with the genuine CP 
violation term, which is not so for the three-neutrino case. We also show that probability 
differences AP a p between CP-conjugate channels have different relative magnitudes in the 
three-neutrino and four-neutrino models. In future, this may provide a way to distinguish 
between these two models. 
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